Discrimination between pure states and mixed states 



Chi Zhang/'B Guoming Wang,-*^'!!] and Mingsheng Ying-'^'0 

^ State Key Laboratory of Intelligent Technology and Systems, 
Department of Computer Science and Technology Tsinghua University, Beijing, China, 100084 

(Dated: February 1, 2008) 

In this paper, we discuss the problem of determining whether a quantum system is in a pure state, 
or in a mixed state. We apply two strategies to settle this problem: the unambiguous discrimination 
and the maximum confidence discrimination. We also proved that the optimal versions of both 
strategies are equivalent. The efficiency of the discrimination is also analyzed. This scheme also 
provides a method to estimate purity of quantum states, and Schmidt numbers of composed systems. 
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I. INTRODUCTION 



In many applications of quantum information, one of 
the important elements which affect the result of quan- 
tum process, is the purity of the quantum states pro- 
duced or utilized. Hence, an interesting and important 
problem in quantum information is to estimate the purity 
of a quantum system [J, |3, S, 0, 0] ■ This problem is also 
strongly related to the estimation of the entanglement of 
multiparty systems [1, 0, Hi ■ 

However, all above references considered this problem 
only in the simplest case of qubits, where both the defini- 
tion and computation about purity are clear. Estimating 
the purity of a general quantum system is still open. In 
this paper, we first consider an extreme situation: given 
some copies of a quantum state, the task of us is to de- 
termine whether the state is pure or mixed. The process 
is called discrimination between pure states and mixed 
states. Then, by counting different results obtained in 
the above discriminations, we offer an effective method 
to estimate the purity of quantum states. The idea of 
discrimination between pure states and mixed states, was 
first mentioned in Ref. ^81]. However, they did not study 
the problem formally and systematically, which is our 
aim in this paper. There are two different strategies to 
design the discrimination: the unambiguous discrimina- 
tion Q , and the maximum confidence discrimination [ld| . 
In the strategy of unambiguous discrimination, one can 
tell whether the quantum system is in a pure state or a 
mixed state without error, but a non-zero probability of 
inconclusive answer is allowed. In this paper, in order to 
simplify the presentation, we use the term "unambigu- 
ous" in a more general sense: it allows the success prob- 
ability to be zero in some situation. On the other hand, 
in the maximum confidence discrimination, an inconclu- 
sive answer is not allowed, and after each discrimination, 
we must give a statement whether the quantum state is 
pure, or mixed. The discrimination is so named, because 
if an answer is given, the probability of obtaining a cor- 
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rect conclusion is maximized. 

It is convenient to introduce some notations here. In 
the Hilbert space iJ®", we use H'^^ to denote its sym- 
metric subspace The orthogonal complement of 
^fyrn Called the asymmetric subspace of iJ**", and 
denoted as Hf^^. We use HH%1^) and ^HfJlJ to 
represent the projectors of these two subspaces respec- 
tively. In this paper, we prove that, given n copies of 
a quantum state p in Hilbert space H, the optimal un- 
ambiguous discrimination and the maximum confidence 
discrimination can be carried out by the same measure- 
ment {Hq = ^H%^^),Ui = The difference 
between these two discriminations comes only from the 
different explanations of the outcomes. In the unambigu- 
ous discrimination, the outcome '0' is an inconclusive an- 
swer, and the outcome '1' indicates that the system is in a 
mixed state. The drawback of the unambiguous discrim- 
ination is that, if the quantum system is in a pure state, 
people always fail to give a confirm answer. However, in 
the maximum confidence discrimination, the outcome '0' 
indicates that the quantum system is considered to be in 
a pure state, and the outcome '1' indicates a mixed state. 

There are two natural assumptions in this paper. 
First, the purity of quantum states is invariant under 
any unitary operation. Suppose the purity of a quan- 
tum state p is represented by p,{p), it must satisfy that 
p{p) — iJ,{UpU^), for any unitary operator U. We also 
assume that, when p is a pure state, p,{p) = 1, other- 
wise < p{p) < 1. For instance, the usually used purity 
of quantum states, p{p) = Tr^(p), clearly satisfies these 
conditions. Second, the priori probability distributions of 
quantum states are also assumed invariant under unitary 
operations. Let us denote the priori probability density 
function as rj(p), then ry(p) = ri{UpU^), for any unitary 
operator U. 

Our present article is organized as follows. In sec- 
tionini we provide the optimal unambiguous discrimi- 
nation between pure states and mixed states. And, in 
section lTVl we provide the maximum confidence discrim- 
ination between pure states and mixed states. We also 
generalize the unambiguous discrimination between pure 
states and mixed states to a "semi-unambiguous" esti- 
mation for ranks of quantum states in section lTTTl which 
can also be seen as an estimation for the Schmidt num- 
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ber of bipartite quantum systems. Finally, in sectionfVl 
we provide a strategy to estimate the purity of quantum 
states. 



II. OPTIMAL UNAMBIGUOUS 
DISCRIMINATION 



and 



Then, 



(7) 



(8) 



In this section, we consider the unambiguous discrim- 
ination between pure states and mixed states. Suppose 
we are given n copies of a quantum state, which is in 
the Hilbert space H . The unambiguous discrimination is 
described by a POVM measurement on the Hilbert space 
i?*^". The measurement is comprised by three positive 
operators, Hp, and II?, satisfying that 



Tr(npp®") = 0, 

for any mixed state 

for any pure state and 



(1) 

(2) 
(3) 



Therefore, if the outcome is 'p', the system is assured to 
be in a pure state; if the outcome is 'm', the system is in 
a mixed state; and outcome '?' denotes an inconclusive 
answer. 

The efficiency of the discrimination is 

P = I Tr(p«"n„0?7(p)dp 

io<M(p)<l 



/ Tr(p®"np)7y(p)dp. 

Ju(p) = l 



(4) 



The optimal unambiguous discrimination is the one with 
the maximum efficiency, and we have the following theo- 
rem. 

Theorem 1 The optimal unambiguous discrimination 
between pure states and mixed states is a POVM mea- 
surement {Hp, Urn, ^?}, such that 



np = o. 



(5) 



where ^(HfJ^) and ^(Hasym) ^''^ the projectors of sym- 
metric subspace and asymmetric subspace of H^^ respec- 
tively. 

Proof. For a mixed state p, whose spectrum decompo- 
sition is p = X]i=i ^i\4'i) {4>i\j and for any n-tuple chosen 
from {1, • • • tt = (tti, • • • ,7r„), where repetition is 

allowed, let us introduce the following two definitions. 



where tt ranges over all n-tuples chosen from {!,••• , m}. 

Because Hp is a positive operator, from Eq.® and 
Eq.©, 



(9) 



for any product state \4'a)- Therefore, when the system 
is in a pure state it satisfies that 



I (gin 



HpIV')®" = 0, 



(10) 



i.e., for any situation, the probability of getting the 'p' 
result is always zero, which means that without loss of 
generality, we can simply let Hp — 0. 

Then, from Eq.([2]), we know that n„i is orthogonal 
to any IV")®", where 1-0) G H. It is known that the span 
space of all jV^}®" is just the symmetric subspace of i?*^", 
which has been denoted as Hf^ [l3]. Thus, the support 
space of n„j must be in the asymmetric subspace H®^^, 
i.e., Ilm < ^{HfJ^j^). The probability of determining a 
mixed state p is 

p{m\p) = Tr(p®"n™) < TT{p®-mfs"ym))- (11) 

Hence, the optimal unambiguous discrimination is the 
measurement {Hp,Hm,H?} given in Eq.((5]). □ 
Under the optimal unambiguous discrimination, when 
the quantum system is in a pure state, the result is sure 
to be inconclusive. If the quantum system is in a mixed 
state p, the probability of receiving an inconclusive an- 
swer is 

p{l\p) = Tr(p®"Hv) 



= y'^pcr(r^). 



(12) 



where is the Gram matrix derived from 

{|0Trj), • • • , |07r„)}, and per(y4) denotes the perma- 
nent of the matrix A, i.e.. 



peviA) = J2Y[A{i,a{i)), 



(13) 



where a ranges over all permutation on n symbols |ll| . 

Let TT is an n-tuple valued in {1, • • ■ , m}. We use 
to denote the number of occurrences of i in tt, where 
« = 1, • • • , m. Because for any two eigenvectors of p with 
non-zero eigenvalues, {(pilipj) = Sij. 



(6) 



F, 



(14) 



i=l 
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where /„t is the n^-dimensional identity matrix. Conse- 
quently, from Eg. p^ 



n! 

7T i — 1 
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SEMI-UNAMBIGUOUS ESTIMATION OF 
SCHMIDT NUMBER 



As we know, the entanglement of a bipartite quan- 
tum system is closely related to the purity of one of its 
subsystems. Whether a subsystem is in a pure state is 
equivalent to whether the total quantum system is in 
a product state. Hence, the measurement given in the 
above section also provides an unambiguous estimation 
for entanglement of bipartite quantum systems. More- 
over, in this section, we will provide a natural generaliza- 
tion, which can be called semi-unambiguous estimation 
of the Schmidt number of bipartite systems. 

The Schmidt number of a bipartite system equals to 
the rank of the quantum state in each of its subsys- 
tems. Hence, estimating the Schmidt number is equiv- 
alent to estimating the rank of quantum states. First, 
let us reconsider the discrimination between pure states 
and mixed states. In the discrimination, the 'm' result 
means that the rank of the state is no less than 2, while 
the inconclusive answer can also be considered as a triv- 
ial conclusion that the rank of the state is no less than 
1. Although the discrimination does not offer the ex- 
act value of the rank of the quantum state, it offers a 
lower bound for the rank. Moreover, the lower bound is 
assured to be correct. In this mean, we can call the dis- 
crimination between pure states and mixed states also a 
"semi-unambiguous" estimation for the rank of quantum 
states. A more general "semi-unambiguous" estimation 
of the rank of quantum states can be defined as a POVM 
measurement on iJ®"- with operators {Hi, H2, ■ • ■ , H™}, 
where m is the dimension of H. The measurement sat- 
isfies that for any quantum state p whose rank is fc, 
Tr(nip®") — 0, for any i > k. Thus, whenever the out- 
come k is observed, we can make sure that the rank of p 
is no less than k. 

Before providing the semi-unambiguous estimation of 
the rank of quantum states. We first introduce some fun- 
damental knowledge about group representation theory 
needed here. For details, please see Ref. [l3| . 

A Young diagram [A] = [Ai, • ■ • , Afc], where ^ Ai = n 
and Ai > A2 > ■ ■ • A/j > 0, is a graphical representation of 
a partition of a natural number n. It consists of n cells, 
arranged in left-justified rows, where the number of cells 
in the ith row is A,-. 



A Young tableau is obtained by placing the numbers 
1, • • • ,n in the n cells of a Young diagram. If the num- 
bers form an increasing sequence along each row and 
each column, the Young tableau is called standard Young 
tableau. For a given Young diagram [A], the number of 
standard Young tableau can be calculated with the hook 
length formula, and denoted by /["^l . In this paper, we use 

tJ'*'' to denoted the rth standard Young tableau, where 
r = l,... 

The Hilbert space i/®", where the dimension of H is 
m, can be decomposed into a set of invariant subspaces 
under operation J7®", for any unitary operation U on H. 
Each of the subspaces corresponds to a standard Young 
tableau tJ'^' , where the number of rows in [A] is no more 
than m. So, we can denote the subspaces as hI^^ , and 
denote its projector as $(i?r'^'). Then, we have 



[A],r 



(16) 



For instance, the symmetric subspace Hfy^ is just one 
of these subspaces, H^^'"^. 

For a quantum state p in H, whose rank is fc, the sup- 
port space of p®" is in the sum of subspaces hI^^ , where 
the number of rows in Young diagram [A] is no greater 
than k. Therefore, a semi-unambiguous estimation of 
the rank of quantum states can be designed as a POVM 
measurement {Hi, • • • , Hm}, such that 



n.- E E*(^^^'') 



(17) 



where /i([A]) is the number of rows in [A]. As said above, 
if the rank of p is fc, Tr(Hip'^") ~ 0, for any i > k. 
Thus, once an 'i' result is observed, we can assert that 
the rank of p is no less than i. For n copies of a bi- 
partite quantum system, through measuring any of its 
subsystems with the measurement given in Eq. (jl7p . we 
can semi-unambiguously estimate the Schmidt number of 
the whole system. 



IV. MAXIMUM CONFIDENCE 
DISCRIMINATION 

In this section, we consider a different strategy 
for determining whether the quantum system is in 
a pure state, which is called "maximum confidence 
discrimination" Td] . 

The discrimination is still a POVM measurement 
{np,H,„}. But when the outcome is 'm', the quantum 
system is believed in a mixed state; otherwise, the out- 
come is 'p', and the quantum state is considered to be 
pure. A maximum confidence strategy is to maximize 
the reliability of the conclusion, i.e., let the following two 
probabilities be maximized. 



p(j)ure\p) = 



//^(p)<i^(^)Tr(p®"n,)dp' 



(18) 
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and 



p{mixed\m) 



/u(p)<i^(p)Tr(p«"n„)dp- 



(19) 



Clearly, Eq. fTSl) and Eq. (fTO|) do not always get maxi- 
mum values at the same time. However, on the assump- 
tions about unitary invariance of ri{p) and /x(p), we can 
prove that there exists a measurement {IlpjIIm} maxi- 
mizing both Eq. p^ and Eg. p^ . as the following theo- 
rem states. 

Theorem 2 The maximum confidence discrimination 
between pure states and mixed states is a POVM mea- 
surement {Hp, Urn} , such that 



where <i>(i?®"j) and ^{Hfp ) are as in Theorem[J\ 



(20) 



Proof. First, we consider the construction of Up. 
From the assumptions that ri{p) ~ rj{U pU^) and p{p) = 
fi{UpW), 

J v{p)TTip^-np)dp 

^(^^^^1^) = .(p)Tr(,-n,)d, 

_ I,ip)=i ^(p)Tr(p^"t/^"np(t/t)g^")rfp 
~ 4(P)<1 ^(p)Tr(p«"f/®"np(C/t)«")dp' 

(21) 

for any unitary operation U. Hence, if Hp maximizes 
Eq.([18l), so does / f7®"Hp(t/1')®"dC/ with respect to the 
normalized invariant measure dU of the unitary group 
U{m). Hence, we can choose the operator Hp to satisfy 
that 



H„ 



C/®"Hp(C/t)®"d[/, 



(22) 



which shows that Hp commutes with any unitary opera- 
tor of the form [/®". Thus, from the representation the- 
ory of classical groups in Ref.[14], Hp can be expressed 
as a linear combination of permutation operators 



Hp = y^QgVg 



(23) 



where G C, cr ranges over all permutations of n ele- 
ments, and Va is the permutation operator derived from 
(T, i.e., 

VMl)\^2) ■ ■ ■ \^n) = \^a^Ma,) ■ ■ ■ I^.J. (24) 

For any state \lp) in the symmetric subspace Hfy^, 
Valf) = \(p), so Ilp\ip) = {J2cr'^'^)\f)^ which indicates 
that 



Up = a^H%l^)(Bll'^, 



(25) 



where H^ is a positive operator whose support space is 



in H^sym, and a = E« 



.acr. Because for any pure state 
p = \tp){-4'\, the support space of p®" is in the symmet- 
ric subspace i?®^, Tr(p®"Hp = for any /x(p) = 1. 
Therefore, the numerator of Eg. pS)) does not change if 
we substitute Hp with a^{Hf^), and the denominator 
diminishes or remains the same. So, the optimal Hp has 
the form of a^{Hf^) for any constant a. 

On the other hand, if we choose $(ff®,y„) as H^, then 
for any pure state p, whose purity p{p) = 1, we have 
Tr(p'^"Hm) = 0, and Eq.([T9|) has the maximum value 
1. To satisfy the condition Hp -t- H,„ = /, let a = 1, 
Hp = $(i7®^). This completes the proof. □ 
It is easy to see that the optimal unambiguous dis- 
crimination and the maximum confidence discrimina- 
tion are the same measurement {Hq = <I>(H^^),Hi = 
^{Hasym)}- The difference between the two discrimi- 
nations is the meaning of the '0' result. In the former 
discrimination, the '0' result means an inconclusive an- 
swer; however, in the latter discrimination, if a '0' result 
is obtained, the quantum system is considered in a pure 
state. 

From Eg. psp . for n copies of a quantum state p, under 
the measurement of {Hq, Hi} given above, the probability 
of receiving a '0' result is 



Po(") 



(26) 



where Ai, • • • , are the eigenvalues of p. As we know, 
the above quantity is the complete symmetric polynomial 
of degree n for {Ai, • ■ • , Am}, which is usually denoted by 
/i„(Ai, • ■ • , A„j). From Ref.[l^, the complete symmetric 
polynomials can be derived from a generating function 



Hm{t) — ^ hk{Xi, • • • , Xm)t^ — 



1 



nr=i(i-^A, 



(27) 



Let A* stand for the maximum eigenvalue of p, then, if 
we have n copies of the states, the probability of judging 
it to be pure can be evaluated as follows: 



Ei=i ni—n i—1 

-\- m — 1 



(28) 



< 



(A*)". 



Then, if the quantum system is in a pure state, Po{n) will 
always be 1, otherwise A* < 1, and poin) will converge 
to zero with exponential convergence rate. 

In sectionHni we discuss the semi-unambiguous esti- 
mation of the rank of quantum states, which is given in 
Ea. p?|) . An open problem is whether this measurement 
also offers a maximum confidence estimation of ranks of 
quantum states, if we consider the result 'i' as a claim 
that the rank of the quantum state is i. 
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V. ESTIMATING PURITY OF STATES 

The maximum confidence discrimination between pure 
states and mixed states provides a natural intuition for 
the purity of a quantum system, i.e., the greater the 
probabihty of getting a '0' resuh, the closer it is to a 
pure state. Hence, by repetitively performing the mea- 
surement, and counting the proportion of '0' results, we 
can estimate the probability of judging the system be- 
ing pure, which, in some sense, reflects some information 
about the purity of the system. However, a more inter- 
esting conclusion is that, no matter how people define 
the purity of quantum states, as long as it satisfies the 
condition of unitary invariant, it can be well estimated 
through a set of maximum confidence discriminations. 

On the assumption of unitary invariant, the purity of 
a quantum state p, ^,{p) = /i(diag(Ai, • • • ,Xm)), where 
Ai,-- - , Am are the eigenvalues of p. Hence, is a 

function of its eigenvalues. Estimating the purity of a 
quantum state p can be reduced to estimating the eigen- 
values of p. The characteristic polynomial of p is a poly- 
nomial, whose roots are the eigenvalues, i.e.. 



det{xl - p) = J[{x - \) 



ajX 



(29) 



3=0 



If we can successfully estimate every coefiicient Oj, j = 
0, • • • , m, the eigenvalues can be estimated by solving the 
equation X^jlo aj^^™ — 0- 

Recall the famous Viete's theorem, it is easy to know 

= eo(Ai, • • • , A„i) = 1 
ai = — ei(Ai, • • • , Am) = — A^ 



for any > 1, if we set er(Ai,-- - ,Am) = 0, when 
r > m. Here, for simplicity, we use e^, hi to denote 
efe(Ai, • • • , Am), hi{Xi,-- ■ , Am) respectively. Then, it is 
not hard to see that 



efe 



hi h2 hz 

1 hi /i2 

1 hi 







hk-i hk 
hk-2 hk-i 
hk-3 hk-2 



hi 
1 



h2 
h2 



(33) 



Clearly, hi = J27Li = Tr(p) = 1. As stated in 
section lTVl for any fc > 2, hk is the probability of receiv- 
ing '0' result, when we measure p®'^ by the measurement 
{Ho = HH%l),Ili = mfslJ}- Therefore, if we have 
N copies of quantum state p, where N is much larger 
than TO, we can estimate the eigenvalues of p in the fol- 
lowing strategy. 

First, separate the N copies into m groups, the kth 
group has kNk copies of the quantum state. Then, oper- 
ate the measurement {Hq = $(i7®^),ni = $(^f®4m)} 
on p®*^ for Nk times in the kth group. Suppose among 
these results, the number of '0' results is Sk, then we can 
estimate po{k), i.e., hk by j^. Then, through Ea. ([55)) . 
we can estimate every e^, wWre 1 < A: < m. Hence, 
from Eq. (|30p . the characteristic polynomial of p, whose 
roots are the eigenvalues we want to estimate is known. 
The task remained for us is to solve the equation given 
in Eq. 



^2 — e2(Ai, • • ■ , Am) — AijAij 

l<2i <22 <n 



ak = i-lfekiXi 



= (-l)™em(Ai 



, Am) = (-!)' 



E - 

l<n<---<ife<n 



, Am) — (— l)™AiA2 • • • A,i 



Here, the polynomial e/c(Ai, ■ • • , A^ 
tary symmetric polynomial of {Ai 
generating function is. 



i) is the m-th elemen- 
, ■ • • , Am} [15], whose 



i=0 



z(Ai 



,Xm)t' ^Yl{l+tX^ 



i=l 



Combined with Eq.(I27|, we have that H{t)E{-t) = 1, so 



fc 

E 



( 1) Gj-hjji—j^ — 0, 



VI. CONCLUSION 



In this paper, we investigate the discrimination be- 
tween pure states and mixed states, which may play an 
important role in further study for estimating the purity 
(30) of quantum states. The discrimination is described by 
POVM measurements {Hq - $(i?f^^),Hi = $(i?®4m)} 
on n copies of the quantum state being discriminated. 
If the '0' result is considered as an inconclusive answer, 
the measurement is the optimal unambiguous discrimi- 
nation. On the other hand, if the '0' result is considered 
as a hint that the quantum system is in a pure state, 
the discrimination is the maximum confidence discrimi- 
nation. We also provide a semi-unambiguous estimation 
for the rank of quantum states, which also can be used to 
estimate the Schmidt number of bipartite quantum sys- 
tems. Finally, we give a strategy to estimate the purity 
(32) of quantum systems. 



(31) 
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